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Abstract 

Properties of the phase space of the standard maps with memory obtained 
from the differential equations with the Riemann-Liouville and Caputo deriva- 
tives arc considered. Properties of the attractors which these fractional dy- 
namical systems demonstrate are different from properties of the regular and 
chaotic attractors of systems without memory: they exist in the asymptotic 
sense, different types of trajectories may lead to the same attracting points, 
trajectories may intersect, and chaotic attractors may overlap. Two maps 
have significant differences in the types of attractors they demonstrate and 
convergence of trajectories to the attracting points and trajectories. Still 
existence of the the most remarkable new type of attractors, "cascade of 
bifurcation type trajectories" , is a common feature of both maps. 
Keywords: Discrete map. Fractional dynamical system, Attractor, Periodic 
trajectory. Map with memory. Stability 
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1. Introduction 

It is commonly accepted tliat fractional differential equations (FDE) play 
an important role in the explanation of many physical phenomena. The phys- 
ical systems that can be described by FDEs, physical fractional dynamical 
systems (FDS), include Hamiltonian systems [Ij; systems of oscillators with 
long range interaction [21 E] ; dielectric ^ and viscoelastic [5J materials; etc. 

The FDSs with time fractional derivatives represent systems with mem- 
ory. Properties of such systems can be significantly different from the prop- 
erties of the systems without memory. As in the case of the regular dynam- 
ical systems, the standard map (SM), or rather the fractional standard map 
(FSM), is a good candidate to start studies of the general properties of the 
FDSs. The first study of the map with memory derived from a differential 
equation, fractional standard map, was done in [H]. References to the prior 
research of the one-dimensional maps with memory, which were not derived 
from the differential equations, can be found in this article. In [6] new types 
of attractors were found for the FSM derived from a differential equation with 
the Riemann-Liouville fractional derivative (FSMRL) and stability analysis 
was performed for the fixed and period two points. 

In this article we present the results of the study of the FSM derived 
from a differential equation with the Caputo fractional derivative (FSMC) 
and compare them with the new and previously [6J obtained results for the 
FSMRL. The results are based on a large, but not exhaustive, number of 
simulations and the continuing investigation may reveal new properties of 
the FSMs. 
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2. Riemann-Liouville and Caputo Fractional Standard Maps 

2.1. Equations and Initial Conditions 
The standard map in tlie form 



Pn+i =Pn- K sin Xn, Xn+1 = x„ + Pn+1 (mod 2n) 



(1) 



can be derived from the differential equation 



oo 



(2) 



n=0 



were e — )■ 0+. 

The equations for Riemann-Liouville and Caputo FSMs were obtained in 
[7] and [8]. The Riemann-Liouville FSM can be derived from the differential 
equation with the Riemann-Liouville fractional derivative describing a kicked 
system 



oo 



oDtx + Ksm{x)J2s{f-in + e)) =0, (1 



< « < 2) 



(3) 



n=0 



were e — )■ 0+, with the initial conditions 



{oDr'x){0+) = Pi, {oDr'x){0+) = b, 



(4) 



where 



oD^x{t) = o/r"x(t) 




(n — 1 < a < n), 



(5) 



D" = d"'/dt"', and qI^^ is a fractional integral. 
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The Caputo FSM can be derived from a similar equation with the Caputo 
fractional derivative 



oo 

CDfx + Ksm(x)J2^{j:-i" + ^))^0, (1< a < 2) (6) 



n=0 

were e — )■ 0+, with the initial conditions 



p(0) = {^Dlxm = {Dlxm = po, x(0) = xo, (7) 

where 



Tin 



(8) 



After integration of equation ([3]) the FSMRL can be written in the form 

Pn+l = Pn - KsinXn, (9) 



a^n+i = T^X^PmK.'l'^-^ + l), (mod27r), (10) 



r a 

^ ^ i=0 



where 



V^{m) = m"-'^ - (m - 1)""*^ (11) 
and momentum p{t) is defined as 

p(t) = oDr'xit). (12) 

Here it is assumed that T = 1 and 1 < a < 2. The condition 6 = is required 
in order to have solutions bounded at t = for a < 2 |6]. In this form the 
FSMRL equations in the limiting case a = 2 coincide with the equations 
for the standard map under the condition xq = 0. For consistency and in 



order to compare corresponding results for all three maps (SM, FSMRL, and 
FSMC) all trajectories considered in this article have the initial condition 
xo = 0. 

Integrating equation ^ with the momentum defined as p = x and as- 
suming T = 1 and 1 < a < 2, one can derive the FSMC in the form 

n-l 

1^^^ V^(n — i + 1) sinxj + sinxn , (mod 27r), (13) 



K 

Pn+l = Pr. 



r(a- 1) 



1=0 



K " 

Xn+i=Xn+Po- ^^:—^^V^{n-i + l)sm.Xi, (mod27r), (14) 



Fa 

^ ' 1=0 



It is important to note that the FSMC ((13), (|14|)) can be considered on 



a torus( X and p mod 27r), a cylinder {x mod 27r), or in an unbounded 



phase space, whereas the FSMRL (([9j), (10)) can be considered only in a 
cylindrical or an unbounded phase space. The FSMRL has no periodicity in 
p and cannot be considered on a torus. This fact is related to the definition 
of the momentum (12) and initial conditions (|4]). The comparison of the 
phase portraits of two FSMs is still possible if we compare the values of the 
X coordinates on the trajectories corresponding to the same values of the 
maps' parameters. 

2.2. Stable Fixed Point 

The SM, the FSMRL, and the FSMC have the same fixed point at (0, 0). 
In the case of the SM this point is stable for K < Kcr = 4. In the case of the 
FSMRL the following system describes the evolution of trajectories near the 
fixed point 

Spn+l = 5pn - K5Xn, (15) 



5 



n+l 



r a 



1 " 



(16) 



2 = 



The system describing the evolution of trajectories near fixed point (0, 0) for 
the FSMC is 

n-l 



5p, 



n+l 



K 



T{a - 1) 



1=0 



K 



5xn+i = Sxn + Spo - — y^V^{n-i + l)SXi. 

1 a' ^-^ 



(17) 



(18) 



i=0 



Direct computations using equations (13), (14) show that the critical curve 
Kcr{a) (see Fig. la) such that fixed point (0,0) is stable for K < Kcr and 
unstable for K > in the case of the FSMC is the same as the critical 
curve obtained from the semi-analytic stability analysis of the FSMRL in [B] . 
This curve can be described by the equation 

ValKcr 



2T{a) 



(19) 



where 



V^i = $^(-l)^+Vi(z). (20) 

1=1 

Taking into account that the equations of the fractional maps (g, ^ and 
(13), (14) and the stability problems (15)-(18) for the maps contain convo- 
lutions, it is reasonable to introduce the generating functions 

oo oo 

X{t) = ^feif and P(t) = X^Pif . (21) 
After the introduction 



i=0 



i=0 



K 



T a. 



5^[(^ + l)"-^-."-^]f, 



(22) 



j=0 



Kit) = + jy^' + - ' ^^^^ 

the stability analysis can be reduced to the analysis of the asymptotic be- 
havior at i = of the derivatives of the analytic functions 

m-'-^^^—, (24) 

P{t)^Po 7 z V (25) 

l-t{l-W^{t)) 

for the FSMRL and 

X{t) = , (26) 



1 



i-t 

for the FSMC. 



tpo + {1 - t)XQ 



Pit) = Po - t ^ -^'^i^W^it) (27) 



2.3. Phase Space for K < K^r and 1 < a < 2 

In what follows, almost all results are conjectures. They were obtained 
by numerical simulations for some values of parameters K and a and then 
verified for some additional values from the corresponding range of the pa- 
rameters' values. 

For the area preserving standard map the stable fixed and periodic points 
are elliptic points - zero Lyapunov exponent. They are surrounded by the 
islands of regular motion. In the case of fractional maps the islands turn 
into basins of attraction associated with the points of attraction or slowly 
diverging attracting trajectories which evolve from the periodic points as a 



decreases from two. For K < K^r and 1 < a < 2 we found no chaotic or 
regular trajectories. Two initially close trajectories that start in the area 
between basins of attraction at first diverge, but then converge to the same 
or different attractors. 

There are significant differences not only between properties of the reg- 
ular and fractional standard maps but also between phase space structures 
of the FSMRL and FSMC There is more than one way to approach an at- 
tracting fixed point of the FSMRL. In Fig. lb there are two trajectories for 
the FSMRL with K = 2 and a = 1.4. The bottom one is a fast converging 
trajectory that starts in the basin of attraction, in which Xn ~ n^^^" and 
Pn ~ (see Fig. Ic). The upper trajectory is an example of the attracting 
slow converging trajectory (ASCT) introduced in [6j in which x„ ~ n""" and 
Pn ~ n^~" (see Fig. Id). The trajectories that converge to the fixed point 
which start outside of the basin of attraction are ASCTs. In the case of the 
FSMC all trajectories converging to the fixed point have the same asymptotic 
behavior: x„ ~ n^^" and p„ ~ n^^"' (see Fig. le). 

In both the FSMRL and FSMC considered on a cylinder the stable fixed 
point (0, 0) is surrounded by a finite basin of attraction, whose width W 
depends on the values of K and a. For example, for K = 3 and a = 1.9 
the width of the basin of attraction is 1.6 < W < 1.7 for the FSMRL ([6]). 
For the FSMC with the same parameters the width is 1.7 < W < 1.8. 
Simulations of thousands of trajectories with pq < 1.6 performed by the 
author, of which only 50 (with 1.6 < po < 1.7) are presented in Fig. If, 
show only converging trajectories, whereas among 50 trajectories with 1.7 < 
Po < 1.8 in Fig 2b there are trajectories converging to the fixed point (0,0) 
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as well as some trajectories converging to the fixed points (0,27r), (0, —2tt), 
and (0, -47r). 

Another significant difference between two fractional standard maps is 
that all periodic points of the SM, except the fixed point (0,0), in the case 
of the FSMRL evolve into periodic attracting slowly diverging trajectories 
(ASDT) (see Figs, la, Ic, le), whereas in the case of the FSMC they evolve 
into the corresponding attracting points (see Figs, lb. Id, If; in Figs. Id, If 
the FSMC is considered on a torus). Presence of only period one structures in 
the phase spaces of two fractional standard maps in the case K = 3, a = 1.9 
(Figs, la, lb) corresponds to the fact that the SM with K = 3 has only one 
central island. Period T = 4 structures in Figs. Ic, Id and period T = 2 and 
T = 3 in Figs, le. If correspond to the phase spaces of the SM with K = 2 
and K = 0.6. Numerical evaluation ([6]) shows that ASDTs which converge 
to trajectories along the p-axis (x — )■ = 0) have the following asymptotic 
behavior: x„ ~ n^^'^ and p„ ~ 

2.4. The FSMRL's Phase Space for K > K^r (I < a < 2) 

For Kcr < K < Kcr2{(^) < 27C, the FSMRL has period T = 2 symmetric 
with respect to the origin stable points with the property 

Pn+l = ^Pn, Xn+l = —Xn, (28) 

which evolve from T = 2 points with the same property of the SM, stable for 
4 < K < 27r. Kcr2 is the upper with respect to K limit of stability of these 
T = 2 points. 



Assumption (x„+i,p„+i) = (-a;„, -p„) = {xupi) in system (d9|, (10)) 
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Figure 1: The FSMRL and FSMC stable fixed point (0, 0): a). The fixed point (0, 0) for 
both the FSMRL and FSMC is stable below the curve K = Kcr{a); b). Two trajectories 
for the FSMRL with K = 2, a = lA, and 10^ iterations on each trajectory. The bottom 
one with po = 0.3 is a fast converging trajectory. The upper trajectory with po = 5.3 
is an example of the FSMRL's ASCT. The value of momentum on this ASCT after 10^ 
iterations is p ~ 0.042; c). Time dependence of the coordinate and momentum for the fast 
converging trajectory from Fig. lb; d). x and p time dependence for the ASCT from the 
Fig. lb; e). x and p time dependence for the FSMC with K = 2, a = 1.4, and po = 0.3; f). 
Evolution of the FSMC trajectories with po = 1.6 + 0.002i, < i < 50 for the case K = 3, 
a = 1.9. The line segments correspond to the nth iteration on the set of trajectories 
with close initial conditions. The evolution of the trajectories with smaller po and for the 
FSMRL with the same K, a, and po < 1-6 is similar. 
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Figure 2: The FSMRL and FSMC phase space for K < Kcr: a). The FSMRL with the 
same values of parameters as in Fig If but po = 1.6 + 0.002i, < i < 50; b). The FSMC 
with the same values of parameters as in Fig If but po = 1.7 + 0.002i, < i < 50; c). 400 
iterations on the FSMRL trajectories with po = 4 + 0.08i, < i < 125 for the case K = 2, 
a = 1.9. Trajectories converging to the fixed point, ASDTs with a; = 0, and period 4 
ASDTs are present; d). 100 iterations on the FSMC trajectories withpo = -3.14-|-0.0314i, 
< i < 200 for the same case as in Fig. 2c (isT = 2, a = 1.9) but considered on a torus. 
In this case all trajectories converge to the fixed point or period four stable attracting 
points; e). 400 iterations on trajectories with po = 2 + 0.04i, < z < 50 for the FSMRL 
case K = 0.6, a = 1.9. Trajectories converging to the fixed point and ASDTs of period 2 
and 3 are present; f). 100 iterations on the FSMC trajectories with po = —3.14 + 0.0314i, 
< z < 200 for the same case as in Fig. 2e {K = 0.6, a = 1.9) considered on a torus. In 
this case all trajectories converge to the fixed point, period two and period three stable 
attracting points. \ \ 



leads asymptotically {n — )■ oo) to 



K 

^1 = T^^^TT^Vaismixi), (29) 
zi [a ) 

Pi = —sm{xi), (30) 



which has been analyzed in The condition of solvability of (29) is 



K > K,r{a) = (31) 

Val 

which means that the stable T = 2 points appear when the fixed point 
becomes unstable. 

An example with K = 4.5 of the evolution of the FSMRL's phase space 
with the decrease of a (from a on the critical curve or from 2 for 4 < K < 2tt) 
is presented in Fig. 3. For large values of a (1.74 < a < 2 for = 4.5) 
the phase space has a pair of stable symmetric T = 2 attracting points 
(Fig. 3a), which in the case 4 < K < 2tt evolve from the centers of the 
corresponding period two islands of stability of the SM. As in the case of 
the fixed point, there are two types of convergence to the attracting points: 
slow with 6xn ~ n~", 6pn ~ and fast with 5x„ ~ n^^^", 6pn ~ n~°. 
For lower values of a (1.67 < a < 1.73 for = 4.5) there appears a couple 
of the non-symmetric stable period two sets of attracting points (Fig. 3b), 
which with the decrease in a (1.63 < a < 1.66 for K = 4.5) transform into 
attracting cascade of bifurcations type trajectories (CBTT) (Figs. 3c, 3d). 
With the further decrease in a the whole phase space of the FSMRL becomes 
chaotic (Fig. 3e). When a is close to one there appears a single chaotic 
attractor, which at the lowest values of a (a < 1.02 for K = 4.5) turns into 
a set of disjoint chaotic attractors (Fig. 3f). 
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Figure 3: The FSMRL's phase space for K = 4.5 > Kcr: a). Period T = 2 stable 
attracting points Xn+i = —Xn, Pn+i = —Pn for a = 1.8 500 iterations on each of 25 
trajectories: po = 0.0001 + 0.08i, < i < 25. Slow and fast converging trajectories, b). 
Two sets of T = 2 stable attracting points Xn+i 7^ —Xn, Pn+i ~Pn for a = 1.71. 500 
iterations on trajectories with the same initial conditions as in Fig. 3a; c). 25000 iterations 
on a single trajectory with a = 1.65, po = 0.3. The trajectory occasionally sticks to one 
of the cascade of bifurcation type trajectories but then always recovers into the chaotic 
sea; d). Time dependence of the coordinate x in Fig. 3c; e). 20000 iterations on a single 
chaotic trajectory with a = 1.45, po = 0.3; f). 7 disjoint chaotic attractors for a = 1.02. 
1000 iterations on each of 20 trajectories: po = 0.0001 + 1.65i, < i < 20. 
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The evolution of the FSMRM's phase space can also be considered for a 
fixed value of a (in this paragraph we consider an example with a = 1.6) 
and K increasing from the value K^r- When the value of K is slightly above 
the critical curve the phase space has a pair of stable symmetric T = 2 
attracting points {Kcr < K < 3.9 for a = 1.6). A couple of non-symmetric 
stable T = 2 sets of attracting points appears for 4 < K < 4.2; CBTTs exist 
for 4.3 < K < 4.4; for K > 4.5 the whole phase space seems chaotic. 

It should be noted that the existence of the CBTTs can't be attributed 
to the fact of the unusual definition of momentum in the derivation of the 
FSMRL. The values of x-coordinates in Fig. 3d represent the solutions of the 
original differential equation ^ independently of the definition of momen- 
tum. 

2.5. The FSMC's Phase Space for K > K^r and 1 < a < 2 

Numerical simulations show that the FSMC also has the period T = 2 
symmetric stable point for K > Kcr- After 10^ iterations the difference in 
the values of xi for two maps is less than 10~^. The order of the difference 
in the values of pi is ten percent. Fig. 4 gives an example of the evolution of 
the FSMC's phase space with the decrease in a for K = 4.5. The intervals of 
the existence of the symmetric T = 2 points, asymmetric sets of the T = 2 
points, and CBTTs are approximately the same for both maps. Trajectories 
converging to the symmetric T = 2 points follow the same law as the trajec- 
tories converging to the fixed point: 5x„ ~ n^~", 6pn ~ n^^°'. It is difficult 
to resolve the CBTTs on the phase portrait of the system (Fig. 4c) and to 
recognize on a zoom but they are very clear on the x of n dependence in 
Fig 4d. The values of the x-coordinates of the CBTTs for the FSMRL and 
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Figure 4: The FSMC's phase space for K = A.b > Kcr: a). Period T = 2 stable 
attracting points Xn+i = —Xn, Pn+i = —Pn for a = 1.8. 1000 iterations on each of 10 
trajectories: po = —3.1415 + 0.628i, < i < 10. b). Two sets of T = 2 stable attracting 
points Xn+i 7^ —Xn, Pn+1 7^ —Pn for a = 1.71. 1000 iterations on trajectories with the 
same initial conditions as in Fig. 4a; c). 30000 iterations on a single trajectory with 
a = 1.65, po = 0.3. The CBTTs can hardly be recognized on the full phase portrait, but 
can be seen on the a; of n dependence in Fig. 4d; d). Time dependence of the coordinate 
X in Fig. 4c; e). 20000 iterations on a single trajectory with a = 1.45, po = 0.3; f). 20000 
iterations on each of two overlapping independent attractors for a = 1.02. The CBTT has 
Po = —1.8855 and the chaotic attractor po = —2.5135. 
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FSMC are very close. For a < 1.6 the phase space becomes chaotic with 
very nonuniform density of points (Fig. 4e). With the further decrease in 
a there appear chaotic attractors and then a set of disjoint and overlapping 
attractors. In Fig. 4f one can see two overlapping attractors for a = 1.02. 
The CBTT on this picture (a dark set) persists even after 300000 iterations 
and overlaps with the chaotic attractor. 

The observed overlapping of attractors and intersection of trajectories 
in phase spaces is a consequence of the fact that the considered fractional 
systems are systems with memory. In such systems the coordinates of the 
next trajectory point are functions of the coordinates of all previously visited 
points on a trajectory. In this coincidence of two points does not lead 

to a coincidence of trajectories, and this property is very different from the 
properties of phase spaces of the regular dynamical systems. 

3. Conclusion 

In this paper we continued the study of fractional maps started in [6]. We 
concentrated on the comparison of two maps at the values of the FSM param- 
eter K < 27C. Even though the simulations were not exhaustive, we were able 
to show that there are significant similarities as well as differences between 
the structures of the phase space of two FSMs. One of the findings is that 
the CBTTs are more common in fractional maps than we originally thought. 
We hope that computer simulations of the equations of the physical systems 
that are described by the fractional differential equations will also produce 
the CBTTs and their origin will have a proper physical interpretation. 
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